Lattice QCD at finite T and \i 

S.D. Katz a * 

a Deutsches Elekronen-Synchrotron DESY, Hamburg, Germany 

Recent results of lattice QCD at finite temperature and density are reviewed. At vanishing density the transition 
temperature, the equation of state and hadron properties are discussed both for the pure gauge theory and for 
dynamical staggered, Wilson and overlap fermions. The second part deals with finite density. There are recent 
results for full QCD at finite temperature and moderate density, while at larger densities QCD-like models are 
studied. 
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1. INTRODUCTION 

Since the first phase transition was seen in lat- 
tice gauge theories in 1981 |1I2| . lattice QCD at 
finite temperature has grown into a large field be- 
coming close to give real, continuum-extrapolated 
predictions for full QCD, which is essential for un- 
derstanding phenomena of the early universe and 
at heavy-ion experiments. As perturbation the- 
ory is only reliable at very high temperatures, lat- 
tice QCD seems to be the only tool to answer such 
questions as the transition temperature between 
hadronic matter and the quark-gluon plasma, the 
equation of state of strongly interacting matter or 
hadron properties at finite temperature and/or 
density. 

1.1. Improved actions 

The formulation of QCD on the lattice is not 
unique. One has different choices for the lat- 
tice action which all give the same continuum 
limit. The simplest choices are the Wilson ac- 
tion in the gauge sector and the Wilson or stag- 
gered action for fermions. Using these actions, 
however may lead to large cutoff effects unless we 
use fine enough lattices. It is often argued, that 
a good compromise is to use improved actions, 
which are computationally (much) more expen- 
sive, but may significantly reduce cutoff effects. 
There are many possibilities for improved actions 
both in the gauge and fermionic sectors. 
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One has to keep in mind, however, that for re- 
liable results, no matter what action is used, sim- 
ulations at least at two different lattice spacings 
(thus different temporal extensions (A^) ), and 
a corresponding continuum extrapolation is re- 
quired. As going to larger Nt is often very costly 
using improved actions, this fact may easily lead 
to the conclusion that using the simplest stan- 
dard actions and a few different A^-s may give 
more precise results than using just one single Nt 
with an improved action. 

1.2. The QCD phase transition 

At high temperature (T) and/or large chemi- 
cal potential (fi) the hadronic matter undergoes 
a transition into a quark-gluon dominated phase. 
At zero fj, this transition happens at T w 270 MeV 
in the quenched case while at T ss 150— 170 MeV 
in full QCD. In the pure gauge case the action 
has an exact Z(3) symmetry which is sponta- 
neously broken in the low temperature phase. 
The Polyakov loop (L) is not invariant under this 
symmetry, thus it is an order parameter of the 
phase transition. In massless dynamical QCD 
we have chiral symmetry, which is expected to 
be broken in the low temperature phase. Corre- 
spondingly, the chiral condensate behaves as an 
order parameter in this case. In the realistic case 
with nonzero quark masses none of these symme- 
try breakings characterize the transition, never- 
theless usually both the Polyakov loop and the 
chiral condensate shows a rapid change (or dis- 
continuity) around the transition. 
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Figure 1. Left: the QCD phase diagram in the m u d,rn s plane The thick solid lines represent second 
order transitions separating cross-over regions from first order ones. Right: the QCD phase diagram in 
the /i, T plane. In the case of two massless flavors, a second order region (dashed line) is connected to a 
first order region via the tricritical point (P). In the realistic case of two light and one heavier flavor, a 
crossover (dotted region) is connected to the first order region via the critical endpoint (E) . 



The transition can be usually located by ob- 
serving the change of these order parameter-like 
quantities. Usually their susceptibilities show a 
nice signal which can be used as a definition of 
the transition point. An alternative possibility 
is to locate the Lee- Yang zeroes of the partition 
function This gives usually more precise val- 
ues for T c . 

An interesting and natural question is, whether 
the transition is of first order, second order or 
maybe only a rapid cross-over. As lattice simu- 
lations are always done on lattices with a finite 
volume (V) where no singularities, thus no real 
phase transitions may occur, the only way to de- 
termine the order of the transition is to examine 
the finite size scaling behavior of certain observ- 
ables, like the susceptibility peaks, Binder cumu- 
lants or the imaginary parts of Lee- Yang zeroes. 

The phase diagram of QCD in the plane of the 
quark masses is shown in Fig. ^ (left). Based 
on universality arguments, one expects a first or- 
der phase transition for three flavors of massless 
quarks and a second order one for two massless 
flavors. In the quenched limit with infinite quark 
masses the transition is of first order as well. At 
nonzero, but finite quark masses the picture is 



more complicated. For intermediate masses we 
expect only a crossover. Whether the physical 
point lies within the first-order or cross-over re- 
gion, is still an open question, however the cross- 
over is much more likely. 

In the latter case an interesting picture emerges 
if we introduce a nonzero chemical potential (/x). 
At zero temperature and large enough /i one ex- 
pects a strong first order phase transition. If 
there is only a cross over at zero /i and finite T 
then there should be a critical endpoint on the 
line connecting the two regions as seen in Fig. ^ 
(right). This endpoint is an unambiguous non- 
perturbative prediction of the QCD Lagrangian 
and has important experimental signatures @|- 

1.3. The equation of state along the line of 
constant physics 

Besides finding the transition temperature 
(which may be not even well defined in the case of 
a cross-over) and/or giving the phase diagram in 
the fi — T plane, one can also give the equation of 
state (EoS) of hadronic/gluonic matter well be- 
low or above the transition T,or \i. All thermo- 
dynamical quantities can be obtained from the 
grand canonical partition function, e.g. for the 
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Figure 2. Left: the continuum extrapolation of T c j ^fa in the pure SU(3) theory using standard and RG- 
improved actions. Right: the pressure and energy density as a function of the temperature for standard 
and RG- improved actions. 



energy density and the pressure: 



computed using the integral method: 



e(T) 



T 2 d(\og Z) 
V dT 



P(T) 



a, 



On lattices with isotropic couplings we can not 
take derivatives with respect to T and V inde- 
pendently, as they are connected via the lattice 
spacing (a) . The combination e — ip can be writ- 
ten using derivatives with respect to a: 



3p 



rp4 



dQogZ) d(logZ) 



da 



(2) 



which can be rewritten in terms of the Plaquette 
(P) and chiral condensate [uu and ss) expecta- 
tion values: 



a& <P> 
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Here all derivatives should be taken along the 
line of constant physics (LCP) where the phys- 
ical values of the quark masses (and not the ami 
lattice masses) are kept constant. In order to de- 
termine both e and p one needs another combina- 
tion, which can be the pressure itself. For large 
homogeneous systems, p oc logZ, which can be 
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dp 



< P > 



da 

- m uT^ <uu> 
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m s — <ss> +/i 
dp 



<9a<9(logdetM) 



dp d(fia) 



(4) 



Here again the integration should follow the LCP 
and the subtraction of T = observables are as- 
sumed. 

In the rest of this review recent results of lattice 
QCD at finite T and/or \x are summarized. In the 
next section quenched results are reviewed. Sect. 
|3l deals with results of dynamical simulations us- 
ing Wilson, staggered or overlap fermions. In 
sect. 0]the recent progress of exploring the /i 7^ 
regions of the QCD phase diagram is discussed. 
Sect. |S] shortly introduces some QCD-like models 
used for studying large chemical potentials and 
finally, sect. concludes. 

2. PURE GAUGE RESULTS 

The simplest case to study is pure SU(3) 
gauge theory without dynamical quarks. Using 
both standard [5] and improved actions |6I7| on 
isotropic and anisotropic [5] lattices the critical 
temperature has been determined in the contin- 
uum limit (see Fig.|2|lcft). The results are usually 
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Figure 3. Left: the chiral extrapolation of T c /y/a for Nf — 2 (squares), 3 (circles) for the p4 and 
Nj — 2 (triangles) for the standard staggered action. Right: the chiral condensate as a function of the 
temperature for several quark masses and extrapolated to m — with Nf — 3 ASQTAD fermions. 



given in units of the string tension and one gets 
T c l^fB w 0.630 - 0.650. The results seem to be 
not completely consistent, however the difference 
comes from possible systematic errors in measur- 
ing the string tension. If one sets the scale using 
the ro parameter P] then a perfect agreement is 
found PHI- The equation of state is shown on 
Fig. |21 (right) for both actions. As the only pa- 
rameter here is the gauge coupling which is used 
to set the scale, the system is kept automatically 
on the LCP. We can see that even at tempera- 
tures as high as 4T C there is still 15% deviation 
from the Stefan-Boltzmann (SB) free gas limit. A 
careful finite size scaling analysis shows that the 
transition is of first order as expected . 

3. DYNAMICAL RESULTS AT n = 

Including dynamical quarks increases the cost 
of simulations by at least two orders of mag- 
nitude. Another difficulty is introduced by the 
quark masses. As direct simulations at the phys- 
ical values of the u,d quark masses are still prac- 
tically impossible, chiral extrapolation is usually 
performed from larger masses. Taking the con- 
tinuum limit is very expensive and it is usually 
not yet done. Therefore we still have to wait for 



final results. 

3.1. Staggered results for T c and the EoS 

The staggered discretization of the fermionic 
action describes four flavors of quarks with equal 
masses. Taking a fractional power of the fermion 
determinant one can also introduce less flavors. 
Using staggered fermions is significantly faster 
than Wilson fermions and another advantage is 
that part of the chiral symmetry is preserved. 
Fig. |3| (left) shows the results of the Bielefeld- 
group on T c using p4 action and 16 3 x 4 lat- 
tices 53- The critical points were defined by 
susceptibility peaks. The critical temperature 
on these N t = 4 lattices in the chiral limit is 
T c « 150(170) MeV for two (three) flavors. The 
MILC collaboration has somewhat different re- 
sults using the ASQTAD action with 2+1 and 3 
flavors [I3j (Fig. right). On 16 3 x 8 lattices 
they see a weak crossover at somewhat larger 
T c w 180 MeV. Note, that although they have 
a larger Nt, the volume they use is much smaller 
than the one used by the Bielefeld group. Thus, 
there may still be significant finite-size effects. 
Clearly, in order to have a reliable prediction for 
T c one has to take the thermodynamical limit and 
carry out a continuum extrapolation. 
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Figure 4. The pressure as a function of T/T c for 
Nf = 3, 2 + 1, 2 (from top to bottom) pA fermions. 



Both of these analyses used degenerate u, d 
quark masses. It has been investigated in |14j 
that breaking the flavor SU(2) symmetry and 
using somewhat different u, d masses does not 
change T c significantly. 

Using standard staggered action the equation 
of state has been determined by the MILC col- 
laboration ^2] on Nt — A and 6. Fig.0|shows the 
pressure obtained by the Bielefeld group with pA 
action on 16 3 x 4 lattices |lfi| . 

Note that both analyses kept am q constant 
which is not an LCP approach as discussed in the 
introduction. In these cases the physical quark 
masses increase as the temperature is increased. 
This leads to a 10-15% suppression of the pres- 
sure at high temperatures compared to the LCP 
approach [T7). 

3.2. Results with Wilson fermions 

There are numerous results for T c using Wilson 
quarks |18I19I2()I21| . Fig. (left) shows the chiral 
extrapolation of T c normalized by the vector bo- 
son mass as obtained by the CP-PACS collabora- 
tion [201 using RG improved gauge- and clover im- 
proved Wilson action on Nt = A lattices. The chi- 
rally extrapolated result is T c = 171(4) MeV, sim- 
ilarly to previously discussed staggered results. 

The equation of state as obtained in |22] is 
shown in Fig. [5] (right) for N t = A and 6 lattices. 



This analysis was carried out along the LCP. 

3.3. First results with overlap fermions 

As the chiral properties of the action may 
be important when we investigate the chiral 
phase transition in the massless limit there would 
be a need to have thermodynamics with chiral 
fermions, i.e. with ones satisfying the Ginsparg- 
Wilson relation [231 - One possible realization is 
the Overlap formalism |24] . Unfortunately due 
to the non-analytic form of the Overlap operator, 
all simulations are computationally extremely ex- 
pensive. In fact, no dynamical simulations have 
been performed before. 

As an exploratory study we implemented the 
standard hybrid Monte-Carlo algorithm for the 
overlap operator |25| . The sign function was 
approximated by the Zolotarev optimal rational 
polynomial approximation 26 . In this approx- 
imation one has to perform simultaneous inver- 
sions which is done with a multi-shift solver |27| . 
As it was expected, the computational costs are 
more than two orders of magnitude higher than 
for dynamical Wilson fermions. 

The first results on a 6 3 x 4 lattice are shown in 
Fig-El The Polyakov loop shows a clear transition 
behavior. This result can of course be considered 
only as a demonstration of the working algorithm, 
for real physical results much larger lattices with 
smaller lattice spacings are required. 

3.4. The order of the transition 

As it has been emphasized in the introduction, 
the order of the phase transition can only be de- 
termined by a careful finite-size scaling analysis. 
In the case of dynamical simulations these anal- 
yses have given little results yet. For Nf — 2 
in the massless limit one expects a second order 
phase transition in the 0(4) universality class. 
However, only simulations with improved Wilson 
quarks could support this expectation |2()j . Stag- 
gered results do not seem completely consistent 
with 0(4) universality |28I29I30| . 

In the case of three degenerate flavors a first or- 
der transition is expected at zero quark masses. 
Thus if we increase the quark masses, we will 
have a critical m c (and correspondingly a crit- 
ical pseudoscalar mass, mps c ) where the first 
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Figure 5. Left: the chiral extrapolation of T c /mv using improved Wilson fermions. Right: the pressure 
for different quark masses on N t = 4 and N t = 6 lattices. 
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Figure 6. The Polyakov loop as a function of the 
(3 coupling using dynamical overlap fermions on 
a 6 3 x 4 lattice. 



order transition ends. This second order end- 
point has been located in standard staggered 
QCD on Nt = 4 lattices jHJ an d the critical 
mass is mps c ~ 290 MeV. However, using im- 
proved p4 fermions, this critical mass goes down 
mpsc ~ 67 MeV [22 indicating that there are still 
large cutoff effects. The picture in the physically 
interesting N j = 2 + 1 case is even more unclear, 
however it is very likely that there is only a rapid 
cross-over 



3.5. Heavy quark free energy and hadron 
properties 

One can study the heavy quark free energy at 
finite temperatures by measuring the correlation 
of Polyakov loops: 

< L X L\ >~ e -VQ*-v\,T)/T (5) 

As an effect of dynamical quarks, string break- 
ing is expected at large separations. This string 
breaking happens at 2-3 times smaller distances 
as we are getting closer to the critical temper- 
ature [34135] . The asymptotic value of the free 
energy at large separations shows a strong tem- 
perature and quark mass dependence. As it is ex- 
pected, for larger temperatures and smaller quark 
masses the asymptotic value gets smaller. 

Measuring hadron masses is not easy at finite 
T, since the temporal extension of the lattice is 
small. One can either use anisotropic lattices or 
measure spatial correlators leading to screening 
masses |3fil37|. 

A recent development is the possibility of 
reconstructing spectral functions from the Eu- 
clidean correlation functions using the Maximal 
Entropy Method (MEM) f53]. For a good de- 
termination, however a large Nt is required. 
Quenched results are already available both for 
anisotropic lattices up to Nt — 54 |2H| and 
isotropic lattices up to Nt = 32 |4U) . 
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Figure 7. The phase diagram of QCD from Nf = 
2 + 1 simulations on N t = 4 lattices. The clotted 
part on the left represents a crossover which is 
separated from the first-order line via the critical 
endpoint. 

4. QCD AT FINITE, BUT SMALL fj, 

Unlike at finite T but zero (i, at finite chemi- 
cal potentials the strategy of lattice simulations is 
not well defined. The reason is that at nonzero \i 
the fermion determinant becomes complex which 
spoils direct importance-sampling based simula- 
tions. 

Recently several methods have been developed 
to extract information for finite p, from sim- 
ulations at zero or purely imaginary /i values 
where the fermion determinant is positive def- 
inite. These techniques are, however still re- 
stricted to finite T and relatively small fi. The 
validity region is approximately p<T. At larger 
chemical potentials, especially at low tempera- 
ture, the only current possibility is to use QCD- 
like effective models which reflect certain proper- 
ties of QCD and which can be solved either ex- 
actly or by numerical techniques. Some of these 
models will be discussed in the next section. 

4.1. Multiparameter reweighting 

One of the possibilities to extract information 
at fi is the Glasgow method |41|. It is based 
on a reweighting in fi. An ensemble is generated 
at fi = and the ratio of the fermion determi- 
nants at finite /i and /i = is taken into account 
as an observable. This method was used to at- 
tempt to locate the phase transition at low T and 



finite fi, however it fails even on lattices as small 
as 4 4 . The reason is the so-called overlap prob- 
lem. The generated configurations (only hadronic 
ones) do not have enough overlap with the con- 
figurations of interest (e.g. in the case of a phase 
transition a mixture of hadronic and quark-gluon 
dominated ones). 

A simple, but powerful generalization of the 
Glasgow method is the overlap improving mul- 
tiparameter reweighting |42|. The partition func- 
tion at finite \x can be rewritten as: 

Z = J VUe~ s ^^ det M(m, \x, U) = 
j VUe- s ^ 0o ' u "> det M(m , [m = 0,U) (6) 



-S g (/3,U)+S g (Po,U) det M(m^ 
det M(mo 



:0,U) / 



where the second line contains a positive definite 
action which can be used to generate the configu- 
rations and the terms in the curly bracket in the 
last line are taken into account as an observable. 
The expectation value of any observable can be 
then written in the form: 

< U >0,m, fl = = 7-5 r (7) 

with w({3, to, (m) being the weights of the configu- 
rations defined by the curly bracket of eqn. (J5J. 

The main difference from the Glasgow method 
is that reweighting is done not only in /j, but 
also in the other parameters of the action (at 
least in 0, but possibly also in m). This way 
the overlap can be improved. If the starting 
point (/3o,TOo,/io — 0) is selected to be at the 
/i = transition point then a much better overlap 
can be obtained with transition points at higher 
\x. One can in general define the best weight 
lines along which - starting from a given point 
(/3o,TOo,/io) ~ the overlap is maximal. This can 
be done e.g. by minimizing the spread of the w 
weights. An alternative way is to diagonalizc the 
covariance matrices for infinitesimal steps |4H| . It 
has been observed that starting from the \i = 
transition point, the above defined best weight 
line coincides with the transition line 44 . 

An interesting question is the goodness of 
reweighting and its dependence on the simulation 
parameters, most importantly on the lattice vol- 
umes. One can define an overlap measure the 
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Figure 8. The phase diagram of QCD from Nf = 
2 simulations on 16 3 x 4 lattices with p4 action 
using Taylor expansion |47|. 

following way: let the overlap measure be 2a if a 
fraction of the configurations with largest weights 
(or their absolute values for complex weights) 
gives 1 — a fraction of the total weight. Clearly for 
true important sampling, when all the weights are 
1, this definition gives 1 for the overlap measure, 
otherwise it is smaller. If we require 2a > 1/2 as 
a condition for acceptable reweighting, then the 
maximal reachable chemical potential, denoted 
by Mi/2 scales with the volume as /U1/2 oc V^ 1 
with 7 = 0.2 - 0.3 from 6 3 x 4, 8 3 x 4, 10 3 x 4 
and 12 3 x 4 lattice simulations [T7]. If 7 < 0.25 
then it has the important consequence that the 
continuum limit can be taken, as the chemical 
potential in lattice units we have to reach scales 
with V' - 25 . 

This reweighting technique made it possible to 
determine the phase diagram on the \i — T plane 
up to fj, « 1.5T in N f = 4 g2] and N f = 2+1 gS] 
staggered QCD on N t — 4 lattices. The transition 
points were located in both cases by finding the 
Lee- Yang zeroes of the partition function. 

In the Nf = 2 + 1 discussed in the in- 

troduction, there is a possibility to have a critical 
endpoint at finite (1e>Te values. Indeed, the fi- 



nite size scaling behavior of the imaginary parts 
of the Lee- Yang zeroes indicated a cross-over at 
fj, = and this turned into a first-order transition 
at larger fi values. After setting the scale with 
T = simulations, the temperature and baryonic 
chemical potential (/ie = 3 • fi) at the endpoint 
can be given in physical units: 

T E = 160±3.5MeV, fi BE = 725±35 MeV. (8) 

The phase diagram and the endpoint is shown 
on Fig. Note that a similar two-parameter 
reweighting was succesfully applied to locate the 
endpoint of the electroweak phase transition even 
on large lattices |4*tj) . 

4.2. Multiparameter reweighting with 
Taylor-expansion 

The use of eqn. © requires the exact calcu- 
lation of determinants on each gauge configura- 
tion which is computationally very expensive. In- 
stead of using the exact formula, one can make a 
Taylor expansion for the determinant ratio in the 
weights |47j (for simplicity assuming no reweight- 
ing in the mass): 



In 



det 



detM(0) J 



E 



(i n <9" lndetM(0) 



(9) 



Taking only the first few terms of the expan- 
sion one gets an approximate reweighting for- 
mula. The advantage of this approximation is 
that the coefficients are derivatives of the fermion 
determinant at fi = 0, which can be well approx- 
imated stochastically. However, due to the ter- 
mination of the series and the errors introduced 
by the stochastic evaluation of the coefficients we 
do not expect this method to work for as large fi 
values as the full technique. Indeed, it has been 
shown in gB] that even for very small lattices (i.e. 
4 4 ) the phase of the determinant is not repro- 
duced by the Taylor expansion for a/i >= 0.2. 

Fig.[S]shows the phase diagram of Nf — 2 QCD 
determined on 16 3 x 4 lattices with pA action using 
this Taylor-expansion technique [47) . The results 
are consistent with the above ones coming from 
the full technique. 
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Figure 9. Left: the phase diagram of QCD from Nf — 4 simulations at imaginary fi with staggered 
fermions [SU]. Right: schematic plot of the phase diagrams (the dotted parts are cross-over while the 
solid parts first order regions) for different quark masses and the line of endpoints (thick solid line). 



As discussed in the previous section Nf = 3 
QCD has a critical quark mass at which there is 
a second order endpoint at \i = 0. If we make 
reweighting also in the quark masses then we can 
keep the system not only on the transition line, 
but also on a critical surface. This way in princi- 
ple the critical endpoint at the physical values of 
the quark masses can be determined |32| . 

4.3. Simulations at imaginary /i 

The fermion determinant is positive definite if 
we use a purely imaginary chemical potential. So 
if the transition line T c (/i) is an analytic func- 
tion then we can determine it for imaginary fi 
values and analytically continue back to real /i-s. 
The analytic continuation is in general impossi- 
ble from just a finite number of points. However, 
taking a Taylor expansion in fi or fi/T one gets: 

(10) 

The coefficients can be determined from imag- 
inary fi simulations. One simply measures T c (/x/) 
for imaginary (ij-s and fits it with a finite order 
polynomial in ni/T. This method has been ap- 
plied for Nf = 2 0H| ,Nf = 4 [SU| and recently for 
N f = 3 |5T| staggered QCD on N t = 4 lattices. 
The Nf — 4 results are shown on Fig. (left). 



The curvatures of the phase diagram in the three 
cases are a 2 =-0.0056, -0.011 and -0.0068. The a 4 
coefficient is very small and it can only be distin- 
guished from zero with very high statistics [51| . 

Similarly to the case of reweighting one can 
study the mass dependence of the fiE endpoint. 
In Nf = 3 QCD for quark masses below the crit- 
ical value this endpoint is at imaginary values of 
H (i.e. negative fi 2 ), which can be investigated. If 
one assumes the same leading order dependence 
on m also for real values of \x (i.e. positive /j 2 ) one 
can in principle determine the location of the crit- 
ical endpoint for the physical values of the quark 
masses (c.f. Fig. El right). The first steps have 
been made in (5Tj and the critical line in the fi—m 
plane has been found for 3 degenerate quarks: 



m e (/i) 
m c (/i = 0) 



= 1 + 



84(36) (A 



(11) 



4.4. The EoS at finite fi 

If we perform reweighting not only along the 
transition line but also below or above it, follow- 
ing the above defined best weight lines then the 
equation of state can be given at finite fj,. For 
calculating the pressure or the energy density one 
simply has to follow the formulae of the introduc- 
tion, but all observables have to be reweighted 
to finite [i according to eqn.J7J). Calculating the 
pressure is done by a line integral in the (3 — fi 
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Figure 10. Left: The pressure difference Ap = p(/i) — p(0) for baryonic chemical potentials mus — 
100,210,330,410,530 MeV (from bottom to top) obtained with N f = 2 + 1 staggered QCD on N t = 4 
lattices. Right: Ap normalized by the Stefan-Boltzmann limit. 



plane. The simplest choice for the integration 
path is to start from deep in the hadronic phase 
(where we get zero contribution after subtract- 
ing the T = contributions) then integrate up to 
some Pi at \i = and then follow a best weight 
line. 

Keeping the system always on the LCP is an 
important aspect of the EoS. However, if we make 
reweighting only in fi and \x then following the 
best weight lines change (3 but not arm, so we au- 
tomatically leave the LCP. There are two possible 
solutions. One could either use also a reweight- 
ing in m^, which makes all computations more 
expensive, or start the reweighting from two dif- 
ferent LCP-s and then keep the system on the 
LCP with an appropriate interpolation. The sec- 
ond technique was used in |52| where the EoS was 
determined for a large range in T and /i. The re- 
sults for the pressure are shown in Fig. ^3 The 
left panel shows the pressure after subtracting 
the fj, = contribution. On the right panel the 
pressure is normalized with the Stefan-Boltzmann 
limit. We can observe an almost universal scal- 
ing behavior. After this normalization the pres- 



sure becomes practically independent of \x in the 
analyzed region. One can get similar results us- 
ing a quasi-particle picture, after setting the free 
parameters of the model using fi = lattice re- 
sults g31. 

The EoS can also be determined using Taylor 
expansion. The pressure can be written as: 




The coefficients can be determined with simula- 
tions at fj, = 0. Fig. ^3 shows the results from 
this method |54j . which are very similar to those 
coming from the full technique. In the vicinity of 
the critical endpoint the nonlinear susceptibility, 
C4 diverges. Thus the convergence radius of the 
p(p,T) function in principle gives the location of 
the endpoint. One can estimate this convergence 
radius by the ratios of the subsequent coefficients. 

The pressure has also been computed with 
the help of ean. (|12|l using quenched configura- 
tions ■ 

Simulations at imaginary \i can also be used to 
determine the a coefficients. A first analysis was 
presented in |5^|. 
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Figure 11. Left: The pressure difference Ap — p{p) — p(Q) using the Taylor expansion method with 
Nf = 2 staggered QCD with pA action on N t = lattices. Right: Ap as a function of /i for different T-s. 



4.5. The factorization method 

Although it has not yet been tested for QCD, 
the factorization method |57| seems very promis- 
ing, so we discuss it in this subsection. The fac- 
torization method was used to study an exactly 
solvable random matrix model defined by: 

Z = J dWe- NT ^ wiw *>detD (13) 

with W being N x N complex matrices and 

"=Ur + „ iW :") (») 

For fi the determinant is complex as in QCD. 
The main idea of the factorization method is to 
use constrained partition functions: 

Z v {x) = j ' dWe~ s \detD\8{x-v) (15) 

with some real-valued observable v. If one per- 
forms simulations at several x values using Z v (x) 
and includes the phase of the determinant as an 
observable then after integrating over x one can 
get < v > for the full theory. The advantage of 
this approach is that each value of the observable 
is sampled with the help of Z v (x) thus solving the 
overlap problem. One disadvantage is of course 
the need of simulations at several x values and 



that for each observable an independent simula- 
tion is needed. The implementation of Z v {x) for 
QCD especially for the most interesting fermionic 
observables (e.g. quark density) is also not trivial. 
However, as this method produced results com- 
pletely consistent with the exact ones for this sim- 
ple random matrix theory j57j . we could expect 
similar success for QCD. The technique is not ex- 
pected to work at large /j, because a reweighting is 
still needed to include the phases of the determi- 
nant and when they become strongly oscillating 
this may make the measurement of the observ- 
ables even with the constrained partition func- 
tions impossible. 

5. QCD-LIKE MODELS 

At low temperature and large densities a rich 
phase structure of QCD is conjectured. Due to 
asymptotic freedom at large densities quarks are 
expected to be weakly coupled. However, as there 
is still a weak attractive interaction they can form 
Cooper-pairs and thus break the SU(3) symmetry 
spontaneously and lead to color superconductiv- 
ity. It is also expected that the ground state will 
be a color-flavor locked (CFL) phase. 

Unfortunately these regions are still unavail- 
able for lattice simulations due to the overlap 
problem and the strong oscillation of the phase 
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Figure 12. The energy-momentum dispersion re- 
lation for the 3+1 dimensional NJL model. The 
energy gap can be seen clearly. 



of the fermion determinant. One can, however 
study QCD-like models, which have positive ac- 
tions and which are expected to reproduce some 
features of high density QCD. In the following 
we briefly discuss the most widely used models to 
study high density phenomena. 

5.1. Two color QCD 

QCD with two colors has a real fermion deter- 
minant, which is also positive definite for an even 
number of flavors even at nonzero chemical po- 
tential. Thus its properties can be studied on the 
whole /it — T plane with standard Monte-Carlo 
methods. Phase transitions have been observed 
both at finite T and [i — and at low T and 
finite 58 52j- An interesting property of two- 
color QCD that it also has an endpoint in the 
fJ> — T plane so it can be used to test any method 
which is used to locate the endpoint of real QCD. 

Unfortunately there is no Fermi surface in this 
model, so it can not be used to study supercon- 
ductivity. 

The modification of hadron properties due to 
the chemical potential has also been investigated 
in this model |6U| . 



5.2. Isospin chemical potential 

If we introduce opposite chemical potentials 
for the u and d quarks then the phases of their 
fermion determinants cancel, so QCD with such 
an isospin fi has a positive definite action. This 
model has been studied in detail recently |61| . 
Interestingly at low densities it behaves very sim- 
ilarly to QCD with real chemical potential, which 
leads to the conclusion that at finite T and low 
density the phase of the fermion determinant does 
not play an important role. On the other hand, at 
low temperatures, just as expected, a phase tran- 
sition happens at fj, w m w /2, as the pions are the 
lowest excitations in this system. This transition 
should not be present in QCD with baryonic \x. 

As in the case of two-color QCD, Fermi phe- 
nomena can not be studied in this model. 

5.3. Four-fermion models 

In principle one could integrate out the gauge 
fields in QCD leading to en effective purely 
fermionic theory. The first non-trivial interaction 
term of such an effective theory is a four-fermion 
interaction. Thus, models with four-fermion in- 
teractions are expected to reflect some features of 
QCD. 

The 2+1 dimensional Gross-Neveu model has 
been studied in |62|. Fermi phenomena have 
been observed, however no gap was found in the 
energy-momentum dispersion relations. The non- 
exponential decay of meson correlation functions 
indicate the existence of massless particle-hole ex- 
citations. 

The NJL model was studied both in 2+1 |53] 
and 3+1 64 dimensions. Interestingly due to the 
different dimensionalities the two models behave 
differently in the large density regions. While in 
both models a qq pairing has been observed, only 
the latter system shows BCS condensation and an 
energy gap in the dispersion relation (see Fig. ll2fl . 
Unfortunately this 3+1 dimensional model is not 
renormalizable, however using a finite cutoff, one 
can set the parameters of the theory to describe 
low-energy observables of QCD. 

5.4. Low energy effective models 

Although the action of QCD at large densi- 
ties is complex, it has been shown in 65 that 
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if one is only interested in the low energy behav- 
ior, i.e. modes close to the Fermi surface, then a 
low-energy effective action can be derived which 
is positive definite (see also jHSl)- Using this low 
energy effective action it has been proven that the 
CFL phase is the ground state of QCD at large 
densities. 

Another approach is to build up an effective 
Hamiltonian based on strong coupling expansion, 
which can then be used to study phenomena even 
in the chiral limit |67j . 

6. CONCLUSIONS 

In this review recent results of lattice QCD at 
finite T and/or \i have been discussed. For the 
pure gauge theory there are already continuum- 
extrapolated predictions both for the phase tran- 
sition temperature and the equation of state. The 
transition is proven to be first order. 

Including dynamical quarks makes all compu- 
tations much more complicated. Therefore there 
are still no chiral- and continuum-extrapolated re- 
sults. However, the transition temperature, the 
equation of state, the heavy quark potential and 
hadron properties are quite well known and all 
necessary, but still missing extrapolations can be 
expected in the near future. 

The case of finite /j, is even more complicated 
due to the complexness of the action. The explo- 
ration of the phase diagram, the critical endpoint 
and the equation of state has just began recently. 
Lots of new results can be expected in the near 
future especially in the small fi region. 

Unfortunately even the latest methods are still 
incapable of solving the complex action problem 
which is most severe for the very interesting large 
fi, small T region. Therefore effective, QCD- 
like models are used to investigate this param- 
eter region. These effective models seem to sup- 
port the assumption on the existence of a color- 
superconducting phase at large densities. 
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